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A THEORY FOR SCATTERING BY DENSITY FLUCTUATIONS 


BASED ON THREE-WAVE INTERACTION 

by 

K. J. Harker and F. W. Crawford 
Institute for Plasma Research 
Stanford University 
Stanford, California 

ABSTRACT 

The theory of scattering by charged particle density fluctuations 
of a plasma is developed for the case or zero magnetic field. The 
source current is derived on the basis of, first, a three-wave interaction 
between the incident and scattered electromagnetic waves and one electro- 
static plasma wave (either Langmuir or ion-acoustic), and second, a 
synchronous interaction between the same two electromagnetic waves and 
the discrete components of the charged particle fluctuations. Previous 
work is generalized by no longer making the assumption that the frequency 
of the electromagnetic waves is large compared to the plasma frequency. 

The general result is then applied to incoherent scatter, and to scatter 
by strongly driven plasma waves. An expansion is carried out for each 
of those cases to determine the lower order corrections to the usual 
high frequency scattering formulas. 
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Introduction 


The scattering of electromagnetic waves by density fluctuations has 
been a topic of general interest for many years. The first derivations, 
given by Booker [l 955 l,and Villars and Weisskopf [ 1955 3 , were based on 
the idea that density fluctuations give rise to dipole-moment density 
fluctuations which in turn cause the familiar far-field electric dipole 
radiation. Most studies since then on scattering use the same basic 
idea. Rosenbluth and Weisskopf [1962] used a technique based on a far- 
field expansion of Maxwell's equations, and a source current consisting of 
a summation over discrete plasma particles. Birmingham et al . [1965], 
although not specifically addressing themselves to the far-field problem, 
showed that this scattering formula must be corrected by a factor equal 
to the refractive index of the scattered wave. 

When the density fluctuations are excited by the random motion of 
charged particles, the scattering is referred to as incoherent scatter. 

The study of incoherent scattering of electromagnetic waves by a plasma 
has been given by a number of authors. Dougherty and Farley [i960], 

Salpeter [i960], and Fejer [i960] independently calculated the cross-section 
for random thermal fluctuations of the electron density. Hagf ors [ 1961 3 
extended the theory to include a static magnetic field. Rosenbluth and 
Rostoker [1962] generalized the theory to take into account departure from 
thermal equilibrium. The subject of scattering by density variations, and in 
particular, incoherent scattering, is thoroughly reviewed by Bekefi [1966]. 

To our knowledge, all of the previous work has been based on the 
high-frequency assumption, i.e. that the incident and scattered electro- 
magnetic waves are much higher in frequency than the plasma frequency. 
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In this paper we generalize this previous work by dispensing with this 
assumption and derive a result which is valid for all frequencies. Of 
course, we still must assume that we are not so close to a resonance that 
we must include multiple scattering effects. 

The source currents responsible for the scattering are determined on 
the basis of two types of interaction, one depending on collective effects 
and one on discrete particle effects. These two effects arise, in turn, 
from the fact that the charged particle distribution function may be resolved 
into two components. One is the spatially averaged part associated with 
plasma waves and collective effects, and the second is the spatially 
rapidly fluctuating component which vanishes when averaged over the 
macroscopic volume. It arises from the discrete motion of the particles 
and is basically a thermal fluctuation phenomenon. 

The mechanism for the collective source current is basically no 
more than a three-wave plasma interaction between the incident and scattered 
electromagnetic waves on one hand, and a scattering electrostatic plasma 
wave on the other hand. The plasma wave may be either a Langmuir or ion- 
acoustic wave. A schematic of the process is shown in Figs. 1 and 2. In 


Fig. 1, the incoming wave mixes with the electrostatic plasma wave 

(<b ,k ) to produce a scattered electromagnetic wave (a> =£h+o k=k+k). 

y 7 v a (3 y } "^7 

In the second version of the process, shown in Fig. 2, the incoming 


electromagnetic wave ,k ) decays into an electrostatic plasma wave 

P 

jiS ) and a scattered electromagnetic wave (tu =co -o> k = k - k 
7 7 v a p y } -XX 


A synchronism diagram showing the dispersion curves of the interacting 


waves and the synchronism parallelogram for the conditions = o)^ ± , 

k = k ± k corresponding to Figs. 1 and 2, respectively, is shown in 

~Q 7 

Figs. 3 and 4 for the case where the electrostatic wave is a Langmuir wave 


and an ion-acoustic wave, respectively 
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The mechanism for the source current arising from discrete particle 


effects is an interaction between the electromagnetic waves again, and 
the synchronous Fourier component of the fluctuating discrete component 
of the electron velocity distribution function. This source current is 
responsible for scattering by unscreened electrons, i . e. scattering which 
does not involve collective effects between the particles. 

Our general mathematical approach is as follows. The far field is 
first determined in terms of an asymptotic expansion of Maxwell’s equation 
(Lighthill, i960). The effects of the two synchronous interactions 
mentioned above are then evaluated by solving the Vlasov equation to second 
order, and using the result to calculate the second order source 
currents. Once the source currents are evaluated, the far field and 
scattered power are determined in terms of products of certain fluctuating 
quantities. If the spectrum of the density fluctuations is known, the 
scattered power is determined by substituting the expressions for these 
products and carrying out the required mathematical manipulations. 

In the case of incoherent scatter, where the density fluctuations are 
not externally driven, but are excited solely by the random motion of the 
plasma particles, it is possible to carry the problem forward to a final 
solution. In this paper we obtain expressions for the product of the 
fluctuating quantities under this assumption, and obtain a closed form 
solution for the incoherent scatter in terms of the unperturbed particle 
velocity distribution functions. 

Since the resulting expression for the scattered power is somewhat 
involved, an expansion in inverse powers of the frequency of the 
incident electromagnetic waves is carried out to gain greater insight into 
the meaning of the results, and to provide a link with the results of previous 


workers . 



The general theory is also applied to the case where the plasma 
waves are so strongly driven by an external source that one can neglect 
the effects of the random motions of the charged particles. Expansions 
are again derived for a high frequency incident electromagnetic wave. 

The results in this paper are based on the assumption that the 
static magnetic field is zero, that the charged particle velocity 
distribution parameters are isotropic in velocity space, and that the 
medium is homogeneous. 
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Theory for Scattering in Terms of Current Sources 


In this section we consider the scattering in general , without 
specifying the current sources responsible. Our system is described 
by Maxwell's equations. 


V X E = - 


Bh 

Tt 


5 e ( i ) ( o) 

V X H = e 0 + + 


(1) 

( 2 ) 


and the Vlasov equation. 


|| + (v * V)f + T] (E + v X B) • || = 0 , (3) 

where JE and JH are the electric and magnetic fields, B is the 
magnetic induction, and are the first and second order current 

densities, f is the electron velocity distribution function, _v is the 
velocity, and T| is the electron charge- to-mass ratio. 

Taking the Fourier transform of Eqs. (l) and (2), and then combining , 
yields the relation 



The first order current is given by 


( 1 ) 

~a 


= mev 

0 ~vc 


' 0 


2 

or e 
ico p ~a 
d a F 


where 

“p = n o e2/me o • 

Equation ( 4 ) then becomes 



e E + —■ — J 

a ~-a di e -a 
a o 



( 4 ) 

( 5 ) 

( 6 ) 

(7) 
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where 


2 ? 

e = 1 - isr /or 

oc pa 


( 8 ) 


Taking the dot product of this equation with respect to gives 

r(2) 


^a‘4 

k • E = — — — - ■ 

-a ~a Vo e a 


( 9 ) 


Substitution in Eq. (7), and solution for E } yields 

OI 




G(C1) j£ ) 

a^a 


4=W “ ’ 


( 10 ) 


where 




( 11 ) 


05 


k (a> ) = — e 1/2 
cc a' c a 


( 12 ) 
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Determination of Far-field Power Flux Density 

We now determine the electric fields in the far-field zone, by 
taking the inverse Fourier transf orm, and then apply essentially an 
asymptotic expansion technique [ Lighthill , i960]. The inverse transform 
of Eq. (lO) is given by 


E(r,t 



G(<d ,k ) 

0 t -jk -r) 

° ^ k 2 -k 2 (cD) 

a a' a 


( 1 $) 


If we replace G by its spatial transform, we obtain 

exp[-jk *(r-r/)] 

E( t ) = — i-r / ■ — p " G(co r^expCjO) t)dk dcD dr' 

(2ny I a £ o e a k 2 -k 2 (cn ) a a ^ a 

' ' f rv rv x rv' 


a cr a' 


Since _r » £ , the integral over k^ can be evaluated in the form 


( 14 ) 


exp[-ik a -(r-r / )] 

k 2 -k 2 (o> ) 

a cr a ' 


2n 2 exp[-ik^(<n^) | r-r' | ] 2rt 2 [-ik Q; (a> Q; )e r . (r-r ' ) ] 


dk 

~a 


\r~£. 


( 15 ) 


where _e^ = _r/r . Finally, the integration over r_ yields 

exp[ j<D t-jk (a) )«r] 

J a J ~cr cr ~- J 


E(x,t) = 



m £ e 
a 0 a 


G[cD k (CD )]dCD . (16) 

~ a 7 —cr cr J a > 


where k (cd ) = k (cd )e . From Eq. (2), the corresponding magnetic field is 

—or a or or —. r 7 

given by 


-B(r,t 



a K, ~oc x cr 


£„(«>) x g[o> ,k (CD )]dCD . 

cr a — a’~a a a 


(17) 
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The time-averaged power flow is given by 


T/2 


P(r)= 1 

T 


im | / dt Re E(r,t) X H*(r,t), 

-♦CO ^ 

-T/2 


( 18 ) 


therefore 


P(r)= lim r f fc 


T// ^ • Ic* exp[j(tu -a/)t-j[k )-k / (o> / )]T} 

dt doi doj' Re; a a a a ~ 


a a 


-T/2 -“ 


r 2 <n (<J3 / )^e _6 e 7 

cr a 0 a a 


(°W a> </) X [~a ( ®y x e (“a^“a>) 


1 


If T is very large, we may take the limit 

T/2 


I: 


xp[j(cu -o/)t dt = 2jt6(a> -o') , 

a a v a cr > 


and Eq. (19) becomes 
P(r)= lim 


e 
T/2 


(19) 


( 20 ) 


X — * 00 


2T( 


L) 3 / 


CO 

dCD 


« rV € e 2 
a 0 a 


* 

g(o> ) x | k (<*> ) X g/o) k (to )H • (2l) 

~\ a'-cr a! \_~a y cr a/~cr cr/J 


Upon substitution of Eqs. (ll) and (l2), and simplification, we obtain 
finally 



An extra factor of 2 has appeared in Eq. (22) because the integration is 
carried out over positive frequencies only. 
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In what follows, it will prove more convenient to write 
Eq„ (22) in its differential form 


SJpfcl = 


30 aaT 


= lim 


e 1//2 cd 2 V 

a a 


a 


7-3— x x 4 2) (v*a' 

(2it )^TVc^e n k ' 




where 


TV -» oo\^ J w '-q a 

0 is the solid angle into which Wave 


a is scattered and V 


the scattering volume. 



Solution of Vlasov Equation 


In the previous section, we derived an expression for the scattered 

( 2 ) 

power as a function of the source current JP ' . In this section we 

a 

determine the latter quantity. In our derivation, we assume an isotropic 
unperturbed electron velocity distribution function, and the absence of a 
static magnetic field. The Fourier transform of Eq. (3) has the form 


df„ 


j (03 -k • v)f +T 1 (E +vXB )-—-H + 

J v a ~a a Wx cr dv 



• 6 ( 



(24) 


where f^ is the unperturbed electron velocity distribution, and the 
subscripts a, 6, and e refer to waves with frequency-wavenumber pairs 
, (tOg,k^), and (o)^,k ) , respectively. Since the incoming wave 
is plane and monochromatic, it has a spectrum of the form 

E = ( 2 *) 4 E 6 (£D -03 ) 6(k-k ), ( 25 ) 


and Eq. (24) becomes 


5f 


5 f 


w J. r — i 

j(VV~ ) V^ ( V~% ) ;‘5T + Z sv. 

S,£ 

where 6,e in the summation run over the values 


(26) 


S = 7 , e = p , 

and 7 refers to the wave for which the synchronism conditions 


0) 

= a> 

+ a> 

a 

P 

7 

k 

~ iL-, 

+ k 

'VC 

“"P 

"7 


hold. 


( 27 ) 

( 28 ) 


( 29 ) 
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We may solve Eq. ( 26 ) iteratively. The first order solution is 
given by 


( 1 ) 

a oj - k • v uae * 

a ~a ~ 


( 30 ) 


where f is the fluctuating part of the solution [ Kadomtsev , 1965] , 
which vanishes when averaged over the macroscopic volume, and satisfies 
the equations 

(<D -k • V ) f _ = 0 , (3l) 

v a ~xx ~~ ua > 

(f (v)f /( v 7 )) = (2k)^ 6(v-v / )6(o -k .v)6(cn -o> /)6(k -k /)f (v) . (32) 

ua ' ua K a —a a a '—a —a ' Qr~' 


Substituting in Eq. ( 26 ), we obtain the second order solution 

f(2) JD V (E + v X B ) • 

a 0) -k . V ZL + ~ 


a -xx ~ 


6,e 


CD -k -v 
e ~e ~ 


+ — 121 — 2 

-k • v c 

a ~a ~ 6 , e 


df 


(^6 + ~ X 


ue 


(33) 
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Second Order Source Currents 


The first order currents, obtained by substitution of Eq. (30) 
into the expression „ 


e 

—a 


/ 4 1} * dv > 


( 34 ) 


need not be considered further, since the contribution from the first 
term on the RHS of Eq. (30) has already been accounted for by Eq. ( 5 ); 
the second term does not contain E as a factor, and therefore does 
not contribute to the scattering. 

Substituting Eq. ( 33 ) into the expression 


-a 


- e /'“ 2) 1 


dv 


gives a second order current 


(33) 


j< 2 > - j< 2) + j , 
—a ~sa —ua 


(36) 


where 


J (2) 

~sa 



d l~e 
' \aT=\ 


£■ • 8 V4 
V* )■ 


(57) 


... v dv (e^+vXB.) 

to -k • v —■ ~ 0 6 

a-~a- 


Bf 

uee 

* “ 57 “ 


(38 ^ 



Source Current from Collective Effects 


We will concentrate first on evaluating jJ^' y the source current 


r(a) 

-sa 


due to collective effects. A partial integration reduces this to 

CO 


,(a) 

-sa 


t ( V V4 2 ~ 


E,+vXB £ 
o ~~ ~~0 . . 

+ - — : > d v 


GD -k • v 

a ~*a ~ 


A second partial integration, followed by expansion of the summation 
according to Eqs. (27) and (28) , yields 

'U 


4a = - 



(E +vXB ) -k 

""6 ~ ^ ~a 

(to -k -v) 2 

v a -a ~ 


, ~7 

to -k -v / 

■y *y (to -k 'v)^ 

1 ' v a ~a 


v + - 


E XB -k 

— y 

(“ 7 ~V^ )( vv^ r 


(39) 


(to) 


(E 'k )(^, + vXB ) E X B 

7 "-a ~3 B y — y ~6 

/ , \2/ \ + (to -k • v)(to -k v} 

(to -k • v ) (to -k -v) v a ~a — ^ -y ~y —J 

a ~a 7 "^7 ■ 7 ' 


k -E 

y y 


(to -k • v ) 


7 “7 L a -a 


(E+vXB )-k E +vXB 

-J & ~ ~ P.„ Zg v + ~P . ~ ~ B 


/ \2 *•' to-k'v 

(to -k -v) a -a ~ 

N ry — ry -*■ ' 


E 

_zg_ 


E • k 
T 


2(k -E ) (k -E ) (k -E )E 

v ~a --g -tot —y ' v + v ~a ~b ~7 


I 


6 (to -k - v ) 2 (to -k «v)(to -k • v)"' 1 (to -k *v)(to -k -v) 2 

p p v a -'•a v p a ~a x p ~p v a "-a 

/ 


In obtaining this equation, we have used the relations 


k *E =0 , 

-a "-a > 


k *e = 0 , 

•^p ®p * 


B = 0 

*7 


(41) 


which follow from the transverse and longitudinal character of the 
linearized electromagnetic and electrostatic waves, respectively. 
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We will find it more convenient to write Eq. (39) in the form 


(2) 

k X k X <r ' 

xx ~xx sa 


Tf eE E k 

■&3-JL- / (v)dv V (v) 
/■' 1 0e w ~s~ 


to 


a 


/ 


( 42 ) 


where 


to lk [(to -k *v)(e *k ) + (k »k )(e *v)] 

v (v) = k x k »-a p - ~rs . -g — 4a ~p =L 

~'~c: rw •~nTV O \ . . O 


~a ~xx ,2 
k 

7 


to (to -k *v) (to -k *v) 
p x a ~<x ~ 7 ~^y ~~ 


(k «k )[(to -k *v)(e .k ) + (k «k )(e *v)] )(k * k )3 V 

/ u \ q --s - ' \ n •'v-x-y/ \«wv • ' -~*-o “v- ' '<w\ ' '*^0 »W' '»-Q — — y^V' «-«< • 


+ 2 


a ^7 p ~p 


p'^a 


N ~na x ~p 


v + 


to (to -k .vP(to -k -v) 
p' a ~*<X v y ~y 


~xx ~p ~-p -a ~p ~a v --p ~7 

p 

(to -k »v) to (to -k *v) 

v a ~a p 7 T 


(k *k )[(to -k *v)e + k (e *v)] [k (e *k ) - e (k *k )] 
a -~y 3 B "-=6 "-Ct ~B ~''P ’**7 


(to -k .v) to (to -k -v) 

3 a •"■a ^ ' p v 7 -’7 


to (to -k *v)(to *k *v) 
p v a ■'a ■ 7 ~7 


(to -k • V V 

' ry -Vy / 


(to -k *v) (e »k ) + (k *k )(e*v) ' v )iU + (j§4 ’v) K, 

P ~P ^ ~P ~xx ~XC ~ p ~p ~ ^ ~P 

. .2 ~ + to (to -k • v) 

VVV** I* 


i§_ 


(k *k ) 2(e *k )(k *k )v 

K ~a -~y' , v ~p ~a v ~a y 


(e.k ) k 
~P ~7 


^P ~P ~ (to -k -v) 2 (to -k • v )^(to -k • v ) (to -k *v) 2 (to -k 'v) 
H H 3 a ~a ~ 3 a ■'■a ~ p ~ 3 a ~a ~ p ~-p ~ 


(43) 


We know from the synchronism conditions [Eqs. (29)] that 

(to -k • v) = (to -k -v) - (to -k ’v) . (44) 

kXkxk=kXkx(k+k)=0. (4p) 

~a ~a -a. ~a — a w p — 7 1 y 

Substituting these into Eq. ( 43 ); an ^ collecting terms, yields 
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r 


CD 




k X k 

~a ~a 


a 


2 

CD k 

P 7 


Tcd -k • v) 
7 ~7 


~ k )3(e*v)k 
"a ^-6 ^7 ~g ~ ~-y 


(cd -k .v)' 

v a ~~a 


k^(e >k ) v (k -k )(k ’k )(e^ *v)v 

cr^ -Hr— W* ~~y ~~Q ~6 ~B ~~ 


(cd -k ‘vY 

K a ~cc 


(co -k *v) c 


(CD -k ’ V ) ^ 

v a -a ~ 


(e 'k )v (k *k )(e *v)v 
V ~"B ~-CC ~ . ~-a B ~B ~ . 

i ~ G 


CD -k *v 
a ~a ~ 


(co -k • v )- 

v a ■'•a 


~P 


(e »v)k 
_ ~B ~~Y 

CD -k *v 

a ~<x ~ 


(k *k )(k -v)e 
K -~g ~y ~ / ~B 

(cd -k -k *v) 

v a ~a ~ P ^ 


(e*k )( k * k )(JL,’ v ) v Cso* k 
f g ^~B ~~y /X,s xx ~-B ~ + B ~~y p ~~ r 


(C6) 


<W*> (“p-Sp-Jt) 


(a) -k • v) -k • v ) 

v a ~xx v p ~13 ■ 



, 4 
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Source Current from Discrete Particle Effects 


The source current due to discrete particle effects is obtained in 
the same manner as the source current from collective effects. We 
expand the summation in Eq . ( 38 ), using only the term corresponding 
to Eq, (27); this is the only term which is dependent on the incoming 
electromagnetic wave, and therefore represents scattering. We obtain 


— ua 


=f~-T.v * d v 

J a ~-a ~ 


V * x V 

A partial integration reduces this to the form 


df. 


uye 


( 2 ) 

j' ' = - iTle If dv 

-ua J / Uye 


E + v X B (E + v X B >k 

~ 3 . Z + - 3 - Z -^! v 


0) -k .V 

a -a ~ 


(tD -k • vV 

v a ~a ' 


(47) 


(48) 


Here again we find it more useful to write this as 


k X k X 

~a ~-a -ua 


j4 e e 

— E 

U) 

a 


/ 


f dv V (v) 
uye u 


where 


(49) 



(^n - ^, * v )e^+!£^( *v) ( U) ^"lL. ,v )( e * k )*(k "k^)(e *v) 

3 ~f> — ~B — g -g ~ —6 ~a A v ~a ~~e 6 ~ 

0) (0) -k »v) + , , ,2 

g v a -a a> (a> -k • 

K g v a -a 


l)‘ 


(50) 
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Scattering Formula : 

We are now in a position to obtain the final scattering formula. 
Substituting Eqs. (42) and (49) into Eq. (2j) gives the equation 


d 2 P 




= lim 


£ 1/2 H 2 e 2 |Ej 2 V 

a J3 


a TV - a, (at^TVc^k 


3rr, r „3^ k ^ 
o a 


1 7)E k 
7 7 


CO 

f t oJ y)dy - 


V S U) 


(51) 


•f 


+ J / f u 7 e ( ^ } 


Since the incoming flux is given by 

and the classical electron radius by 

r 0 = e 2 /(4jt e Q mc 2 ) , 


(52) 


(53) 


the scattering formula can simply be written as 


S 2 p 


= lim 


iv « 1/5 


ST&ir " ““ ^TJ4 v . o 

a TV-, „ k a P' 


11 E k lf n (v)dv V (v) + j/f dv V (v)( 2 

'v -v I ,bw s - J / n<v<=» — ^ir^ 7 1 


(rj 


uye m 


(54) 


Expanding the squared term gives the equation 


d 2 p 


a 


= lim 
TV 00 


a 


/e \ 1/2 1 

f 


’i 2 k ^ c 5^ ] ;i/ f o e ( ^ s w* 

k *—co 


2 Re jT1k 7 / f 0e ( _v)d_v V g ( 


*>) - ( / [v w * E r ] %&*) 


(55) 


+ 

—CO 


Co 

/ d.vdv'Tf (v) f* (v' )] V (v) •V*(v / ) ( . 
J |_ uye — u/e ^ J —u— ~u w ' ( * 
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This equation is the general scattering formula we have sought to 

derive. If one knows the spectrum corresponding to Ie 1^, f* E ' and 

7 7 uye 7 5 

If* I ^ and the unperturbed velocity distributions then the scattered 
U7e 

power is determined. 
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Incoherent Scatter 


We will now take up the case of incoherent scatter ^ where it is possible 
to evaluate Eq. ( 55 ) explicitly. In this case the assumption that the 
charged particle motions are random allows one to evaluate the products 
of the fluctuating quantities in the equation. In the appendix we show 
that these products are given by 
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Substituting these expressions, and integrating over v/, reduces Eq. (55) 
to the form 
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Because of the delta function and Eq. (44), we can replace the 
factor (01 -k -v) by(oi -k *v 
this out, along with the application of Eq. (45), yields the simpler equations 


)in the definitions for 


^2 


and L_ 
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Equation (59), along with the definition given by Eqs. (46), (60), ( 63 - 65 ), 
is our final result describing scattering by random density fluctuations. 
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High Frequency Expansion for Incoherent Scatter 

Equation (59) is a final result in the sense that it specifies 
completely the scattered power once the unperturbed velocity distribution 
functions are known. It will be useful , however, to expand this equation in 
powers of 1/03 , in order to interpret the meaning of the result, and to 

P 

compare with previous work. 


0 ) 


Let us first expand L in Eq . (60) to second order in 

~'l P 

yields, upon application of Eq. (4-4), the relation 
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If we substitute into this equation the vector identity 


v = — 7 T [k (k -v) - k X k X v] , 
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(67) 


and collect terms, we obtain 
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+ terms in (k X k X v) . 
Finally, substituting into Eq . (60) gives the result 


( 68 ) 
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is the charged particle susceptibility. In obtaining this result we have 


ignored the terms in V (v_) containing (k X k X v} , since these give 

~g y -~-y 

rise to terms of order (v^/c) ,and higher, when dealing with isotropic 
electron velocity distribution functions. 

Let us now expand in Eq. (63). To second order in 0 ) ^ we obtain 

/ oj ( r fi k ~ 

dv f_ (v)6(u> -k *v) — k X k X '< e_ + (e *k W-/ e • vVk 1 — — + ' 1 

*“ 0e v ~' 7 ~y 03 ^ -o' ~a j~-p ~ar \ p —j~y J ^2 


(k - k )(e *v)v 

'-a - J e / v ^p ~ ~ 


We again expand jv according to Eq. (67)>a nd ignore the component perpendi- 
cular to _k^ . This yields, after replacing (jc^ • v) by , the result 
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A similar procedure for JL yields the expression 
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We now substitute the expressions for and L' into Eq. (59) 

to obtain 
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Using the identities 
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and rearranging, gives the final form for the expanded incoherent 
scattering formula as 
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When the frequency of the incident electromagnetic wave, 

0)^ } becomes very large compared to and cl>^ , the high frequency inco- 

herent scattering formula [ Bekef i , 19 66] , 
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is retrieved. In the case of backscatter, where ( e • k^) = 0 } Eq. (77) 
reduces to the even simpler form 
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Scattering in Case of Strongly Driven Plasma Waves 

In the case where the coherent waves are so strongly driven by an 

external source that the random fluctuations of the charged particles can 

be ignored , we may ignore the terms involving f in Eq. ( 55 ), and 

uy e 

the scattering is then given simply by 
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As in the case of incoherent scatter, it is useful to determine the 
behavior for a high frequency incident wave. Substituting Eqs. (60) and 
(69) into this formula reduces it to the form 
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we 


obtain the standard high frequency formula [ Bekef i , 19663 given by 
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The ( e c / € p^ 1// ^ correction to Eq. ( 83 ) y contained in Eq. (80), was Previously 
derived by Birmingham et al . [ 1965 I by a different method. 
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Summary 


A general theory for scattering of electromagnetic waves by density 
fluctuations in a plasma has been carried out. The general scattering 
formula is given by Eqs. (55), (46), and (50) . Its application to incoherent 
scatter is given by Eqs. (59) , (46), (60), (63-65), and to scatter by 
strongly driven plasma waves by Eqs. (46) and (80). The theory generalizes 
previous high-frequency theories in that it is valid for all frequencies 
of the incident and scattered electromagnetic waves. It does assume, 
however, a zero magnetic field, isotropic unperturbed charged particle 
velocity distribution functions, and the absence of multiple scattering. 

An expansion for both incoherent scatter and scattering by strongly 
driven plasma waves in inverse powers of the frequency, cd 1 of the 

P 

incident electromagnetic wave has been carried out [Eqs. (77) and (8l) 
respectively]. These expansions show that two types of lower order 
corrections must be applied to the high frequency theory as the incident 

electromagnetic wave frequency approaches the plasma frequency. The 

2 

first type of correction is of order (cd /CD ) and must be applied 

P p 

irrespective of the value of the difference frequency, CD } between 
the electromagnetic waves. The second correction is of order (cd^/cd^) , 

and is clearly of importance only for scattering by the Langmuir waves. 

These lower order corrections disappear for the case of backscatter. 

As the frequency, cd^ , of the electromagnetic wave comes closer to 
<d^ , then of course it is necessary to use the full theory [Eqs. (59 ) , 

(46), (60), and (63-65), or Eqs. ('46) and (80)]. It is important to note 
that the full theory has non-vanishing higher order corrections for the 
backscatter case, even though the lower order corrections mentioned in 
the previous paragraph disappear. 
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Appendix 


In this section will be derived the expressions for the space-time 
averages given in Eqs. (56-58). Our first step is to derive a relation 
for the averages in terms of the ensemble average of the Fourier components. 
According to Parseval's theorem, the average of the product of two variables, 
A(jr, t) and b(j:, t), is given by 
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The ensemble average, on the other hand, is given by 
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All of the cases studied in this paper have the property that 
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Equating the two averages given by Eqs. (A.l) and (A.k) yields the desired 
relation 
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Substituting this into Eq. (A. 5 , ),with A and B equal to f u ('v) and f (v*), 
respectively^ yields Eq . ( 56 ) immediately. 

In order to prove Eqs. (57) and (58), we will need to use the 
linearized Poisson equation 
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If we multiply this equation by f* /(v 7 ), take the ensemble average 

U 7 e ' ' 

and assume that the ion and electron motions are uncorrelated, we obtain 
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Substituting Eq. (32) shows that 
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Substituting this result into Eq. (A. 5 ), with A and B equal to E and 
* 

f , respectively, yields .Eq. (57). 

If we multiply Eq. ( A. 7 ) by its complex conjugate, and take the 
ensemble average,we obtain 
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Substituting Eq. ( 32 ) allows us to determine C (EE : ). If this is in turn 


substituted into Eq . (A. 5 ) } with A and B equal to E, we obtain Eq. ( 58 ). 
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FIG. 1. Mixing of an incoming transverse wave (k ) 
and an electrostatic wave (k^) to produce a 
scattered transverse wave (k^) . 



into an electrostatic wave (k ) and a scattered 
transverse wave fk 



FIG. 3. Synchronism diagram for the interaction of two transverse 


waves and a Langmuir wave. 
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FIG. 4. Synchronism diagram for the interaction of two transverse 
waves and an ion-acoustic wave. 
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